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ABSTRACT 

Second-order steady wave forces on a restrained, submerged, 
circular cylinder are studied experimentally as a function of 
cylinder depth and wave frequency. The experimental results are 
compared to forces calculated from a linearized potential theory. 
It is found that in all cases the experimentally observed forces 
are higher - in some cases by an order of magnitude - and are not 
attenuated with depth as rapidly as predicted by the theory. 
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1. INTRODUCTION 



Since very early in written history warriors have 
attempted to escape the sea's surface and utilize submarine 
warfare. Herodotus (460 B.C.), Aristotle (332 B.C.), and 
Pliny, the elder (7? A.D.) mention determined attempts to 
build submersibles . Alexander the Great (35^ to 323 B.C.) 
is the first person known to have descended into the sea in 
a vessel of any kind.^ 

Since that time submarines have been invented which 
were propelled by oars , sails , treadles , hand-operated screws , 
clockwork, springs, steam stored in tubes, chemical engines, 
compressed air, stored gases, electric motors, and nuclear 
power. Over three hundred years ago Nother Shipton, a famous 
English prophetess , predicted the coming of the submarine when 
writing, '‘underwater man shall walk, shall ride, shall sleep, 
shall talk." Her prophesy has come true, but in carrying 
out her prophesy man has discovered that even though he may 
be an underwater man in many senses, he has not completely 
escaped the forces of surface waves. 

Man was early to recognize the forces of the sea to be 
powerful, at times destructive, and inevitable. He set out to 
study them, attempt to predict them, and to design vessels 
which would withstand and take advantage of them. Studies of 
forces on vessels penetrating a free surface are too numerous 



to mention here. 



The United States Navy commissioned the USS Holland, its 
first submarine on April 11, 1900, The military importance of 
submarines has increased at a staggering rate since that time. 
More recently, commercial and scientific interests have awakened 
to the vast resources under the sea, hence the use of non- 
military submersibles has begun to increase at a rate that may 
exceed that of its military cousin. As more and more use is 
piade of undersea vessels , they become larger — the largest 
known maneuverable undersea vessel to date is about ten meters 
in diameter, greater than one hundred meters in length, and 
displaces in excess of six thousand tons — and yet they must 
be precisely controlled. 



In order to be able to provide this precise control, the 

forces acting on these vessels must be known. One source of 

such forces is surface waves. A number of studies have been 

made of the wave induced forces on underwater bodies, but most 

of these ‘studies have been concerned with an underwater body 

moving with some translational velocity with respect to the 

2 

fluid. Examples of these studies are some of Havelock's 

3 4 

papers, Kim's J paper, and Salvesen's paper. 



The present work will concern itself with a body that is 
not moving, but which is rigidly restrained. Kulin^ studied the 
wave forces on submerged cylinders and plates by considering 
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the effect of a single intumescence which subjects the fluid to 
acceleration followed by deceleration. He studied (l) the 
extent local acceleration causes the drag coefficient to depart 
from steady state values, (2) the extent that friction affects 
potential flow values of inertial coefficient, and (3) the extent 
of force dependence upon preceding flow history, as well as on 
instantaneous conditions. He found that the forces could not be 
represented by a single universal total resistance coefficient, 
and that the ratio of the distance of fluid motion to object 
size was a significant parameter related to vortex formation 
during fluid deceleration. 

In 1962, Ogilvie^ applied and extended the form of 

7 

solution that had been developed by Ursell in 1950 to calculate 
the first-order oscillatory force and the second-order steady 
force in the following situations: (a) a submerged cylinder is 
restrained from moving under the effect of incident sinusoidal 
waves; (b) a submerged cylinder is forced to oscillate sinusoidally, 
in otherwise calm water; (c) a cylinder, which is neutrally 
buoyant, is allowed to respond to first-order oscillatory forces . 

He proves that subject to his assumptions a knowledge of the first- 
order potential supplies information sufficient to solve these 
problems. This paper deals with an attempt to experimentally 
verify his results for case (a) , the restained cylinder. 



Assume that a circular cylinder is located under a free surface, 



(see figure 1.) such that the surface of the cylinder is 
specified by 



S(x,y) = x 2 + (y + hf -a 2 =0„ (l) 

The undisturbed free surface is taken as the x-axis and 



the instantaneous free surface will be specified by 

y - Y(x,t) = 0 o (2) 

The required velocity potential, 23(x,y,t), satisfies 

= ° for y<Y(x,t), (3) 

and S(x,y) > 0,' 

® a Y x - ffi Y + Y 4 = 0 on y = Y(x,t) , (4) 

gY + E t + |(2£ + 23*) = 0 on y = Y(x,t), (5) 

and 

iD x S x + 23 y S r + S f = 0 on S(x,y) = 0, (6) 



plus appropriate conditions at infinity. The time average, Z u, (t) , 
of the second-order steady force can be shown to be 

z*Tty*= x'Ttyh i? E rty t 

it 1 

= e iQ (m y ) x J ^^d© (?) 

-ir 

or if we let 2!(x,y,t) be the real part of a function of 
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a complex variable, f(z,t), where z = x + iy, then 

*jr r 



&Ktj - -}±afj e lQ f '(z,t)f ' (z.t)^^© 






~2kh 

= -i27Tpg(h T 0) e Q(ka,kh). 



k is the wave number defined as 



k =0 J /g. 



and HO is the amplitude of the incident waves. 



Define steady vertical force coefficient,^, as 



e = ?rtr • 

277js>g(HO) 

From (8') and (10), it is seen that 



fi- 



-2kh 



Q(ka,kh) , 



It is this quantity obtained by Ogilvie that will be 
compared -with experimental results in Section IV. 



( 8 ) 

( 8 *) 



(9) 



( 10 ) 



( 11 ) 
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Figure 1, Geometry of the problem 
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II. TEST APPARATUS AND PROCEDURE 



II. A. TEST APPARATUS DESCRIPTION 

This experiment was conducted in the wave tank in the M.I.T, 
Department of Naval Architecture and Karine Engineering Hydro- 
dynamics Laboratory. The working section of the tank is ten 
feet long, one foot wide, and six inches deep. The wave generator 
is a vertically oriented pusher, hinged at the bottom, and driven 
by a push-rod connected off center on a f lywheel , which is in turn 
driven by a variable speed electric motor. The observation section 
was located about equidistant from the two ends. An energy 
absorbing beach made of rubberized horse hair was located at 
the end opposite the wave generator. This beach was not completely 
successful in preventing reflected waves. This will be considered 
in Section II. 3. 

Wave height was measured by a capacj tance probe mounted on 
a carriage which was driven by a variable speed electric motor 
and ran on tracks along the top of the tank. The output of the 
probe was fed through a buffer amplifier to a Hewlett Packard 
Model 7702 B recorder. The carriage was moved at a constant 
speed before each force observation to obtain data for cal- 
culating the reflected to direct wave ratio and then placed 
at a predesignated point as close as possible to the cylinder 
to measure wave height during the period of force measurement. 
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The cylinder itself, a one inch diameter lucite rod 11 15/l6 
inches long, was suspended below the surface by a 3/32 inch 
diameter rod tapped into the cylinder at points l/8 inch from 
each cylinder end. Any effect on vertical forces by these thin 
supporting rods was assumed to be negligible. The l/32 inch 
clearance between each end of the rod and the tank wall is 
considered to be smaller than any wall boundary layer thickness , 
and it is felt that there was no appreciable departure from the 
two dimensional condition. The tops of the supporting rods 
were threaded and attached to a yoke which was restrained so 
that its only degree of freedom was in the vertical direction. 

Vertical force on the yoke was measured by two Schaevitz 
Engineering Kodel Ko. FTA 3 Force Transducers, one on each end. 
The electrical outputs of these transducers were connected and 
processed as shown schematically in figure 2, As seen from the 
figure the output voltage was directly proportional to the time 
integral of the steady force component, a. 
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Figure 2, Functional force measurement circuit 



II. B. PROCEDURE 



As mentioned in Section II. A., the energy absorbing beach 
was not completely successful in preventing reflected waves and 
it was considered necessary to take these reflected waves into 
account. A method of doing this by the use of a standing wave 
ratio, p, is developed in Appendix A. That result is given in 
(A7) and is- repeated here. 



The values of HX and HE used in (A8) were obtained by 
allowing the capacitance probe carriage to move at constant 
speed along the tracks as mentioned in Section II. A. and observ- 
ing the maximum and minimum values of wave height recorded. 

The wave height recorder was calibrated by changing the 
still water level by a known amount and observing the change 
in recorder level, Eecause of the sensitivity of a capacitance 
probe , the calibration was repeated frequently and the appropriate 
calibration factor was used in each calculation of B . 



Force measurement calibration was accomplished by first 
opening the shorting switch shown bypassing the integrating 
amplifier in figure 2 . for a known period of time , TD , and 
observing the integrated drift voltage, ED. The switch was then 
closed and a known w^^t. ^ t,toq •nlar'orl nn bVin nvl i nd rt* sunnortinc 

• ■ w — O'* f f — — — -*/ XX “W 




( 12 ) 



- 14 - 



yoke 0 The switch was reopened for a known period of time, TS , 
and the integrated signal voltage, ES , was observed. The cali- 
bration factor, AE, was then determined by the folloxring formula; 



AE = Y/E 



( 13 ) 



where 



E = (ES/TS) - (ED/TD) . 



(1*0 



This factor was found to remain constant at 0,036 pounds per volt. 

The unknown steady wave force, l w, (t) , was found by the 
following equation , using values of ES , TS , ED , TD , HM , and HX 
observed for various cylinder depths, h, and incident wave 
frequencies, F, 



y#0was calculated using (10) and ( 15 ) a^id was calculated by (12), 

Calculations were performed at the M.I.T. Computer Center 
using the program in Appendix B. 



T<ro*= (E) (AE) 



( 15 ) 
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III. RESULTS 



Experimental data is tabulated in Appendix C„ Values 
of ft calculated from this data using (12) are tabulated on 
the following page. 

These values of are presented graphically in figures 
3. through 9. along with curves for ka = 0.2 and 0.5 from 
Ogilvie's theory. 

In drawing the experimental 0 curves , curves for the next 
higher and lower values of ka were considered as well as /3 
values for the curve being drawn. As a result the experimental 
curves for some ka values do not pass through all data points. 

The individual /3 curves from figures 3« through 9* are 
reproduced in figure 10. for ease of comparison. 
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Figure 3. Steady vertical force coefficient, ka = 0.25 
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Figure 4. steady vertical force coefficient, ka 
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Figure 5. Steady vertical force coefficient, ka 
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Figure 6. Steady vertical force coefficient, ka =0.4 
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Figure 7. Steady vertical force coefficient, ka = 0.^5 
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Figure 8. Steady vertical force coefficient, ka = 0.5 
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Figure 9* Steady vertical force coefficient, ka = 0.55 
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Figure 10. Steady vertical force coefficients 
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. IV. DISCUSSION OF RESULTS 



Reference to figures 3» through 9» shows that the 
experimental curves of steady vertical force coefficient are 
in all cases higher — in some cases more than an order of 
magnitude — than the theoretical curves , and in all cases 
except for ka = 0.25 and 0.3 the slope is less steep than the 
theoretical curves. 

It is also noted that the experimental curves all have 
inflection points between the third and fourth data points — 
at an h/a ratio of about 4, From figure 10. it is easily 
seen that this inflection is much more pronounced in the curves 
for smaller values of ka — longer wave lengths. Since for the 
longest wavelength investigated the tank depth was over a half 
wave length it was expected that it would have an infinite 
effective depth; however, a finite depth effect seems to be a 
plausible explanation for these curve inflections. 
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V. CONCLUSIONS AND RECOlliENDATIONS 



From Sections III and IV it may be concluded that Ogilvie's 
potential theory does predict the general behavior of the steady 
vertical force coefficient, but that in all cases it predicts 
a force smaller than that measured experimentally and that the 
forces will be attenuated, more rapidly with depth than that 
experimentally observed. 

Because of the curve inflections noted in Section IV and 
since smaller values of ka are more likely to be encountered at 
sea, i.e. , a one hundred meter wave length and a ten meter 
diameter cylinder give a ka of 0 . 31 ^ — most waves of interest 
are longer and most diameters are smaller, hence most ka's of 
interest are smaller, it is felt that an extension of this 
experiment .to smaller values of ka in facility with an infinite 
effective depth and no reflected waves would be rewarding . 
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APPENDIX A 



Theory Modification to Account for Reflected Waves 

g 

Using the 'circle theorem' of Milne-Thompson, the complex 
potential function, f(z,t), in (8) can be represented by 
4(<rt + kz) -kh ijcrt + ka* 



f, (z,t) = A |e + e el (z + ih), 

where 

A = (HO) g/<r. 

Using (Al) , (8) and (10) , is found to be 
-2kh 

ft = e (ka) I, (2ka) . 

The complex potential function for a reflected wave 
traveling in the opposite direction with an unknown phase 
angle, H 7 , and a magnitude of p times the direct wave can be 
given by 



(Al) 



(A2) 



(A3) 



iT 


-i(crt - kz) -kh i 


1 

A 


f a (z,t) = pAe 


e + e e 


(z - ih) 

’ 'J 



(A4) 



Now the complex potential function for a direct and reflected 
wave existing simultaneously is given by 

f(z,t) = f, (z,t) + f z (z,t), (A 5) 

Using (A5) , (8) and (lO),ySO is found to be 

-2kh 



/SO = (1 + p 2 ) e (ka) I, (2ka). 



(A6) 



Consequently 



0 =P°/( i + p 2 ) 



(A?) 



where ft would be the steady vertical force coefficient if 
there were no reflected wave present , 00 is the steady vertical 
force coefficient with a reflected wave present, p is the 
standing wave ratio given by 



and HX and HM are the maximum and minimum wave amplitudes 
that exist as a result of the summation of the direct and 
reflected waves. HO used in (A2) is the average wave amplitude 
given by 



The complex potential functions given by (Al) and (A4) do 
not satisfy the free surface conditions (4) and (5) » hence 
(A3) and (A6) can be expected to give reasonable results only 
if ka « 1 and/or 2kh » 1. However, it is assumed that (A7) is 
valid throughout the interest range of ka and kh values. 



P = HX - HM 
HX + HM 



(A8) 



HO = HX + HM 
2 



(A9) 
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APPENDIX B 



Computer Program for Steady Vertical Force Coefficient Calculation 
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APPENDIX C 



Experimental Data Tabulation 
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